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Abstract 

It is shown that certain classes of differential equations arising from the mod- 
elling of chemical reaction networks have the following property: the state space 
is foliated by invariant subspaces each of which contains a unique equilibrium 
which, in turn, attracts all initial conditions on the associated subspace. 
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1. Introduction 

There are difficult and interesting open questions about allowed asymptotic 
behaviour in systems of differential equations arising in the modelling of chem- 
ical reaction networks (CRNs for short). The main goal in this area is to make 
claims about the behaviour of these systems which are as far as possible in- 
dependent of the particular choices of functions or parameters which describe 
the rates of reaction or "kinetics" . Classical results in this direction [l|, 0| rely 
strongly on the choice of "mass action kinetics" leading to particular polynomial 
differential equations. Mathematically, such results involve proving that solu- 
tions of certain parameterised families of polynomial differential equations have 
certain asymptotic behaviours regardless of the values of the parameters, but 
provided these have fixed sign. However when the kinetics is constrained only 
by loose qualitative laws, the family of possible differential equations describing 
a reaction network becomes much larger, and results become fewer. Here, we 
provide a general result based on the theory of monotone dynamical systems 
[H, HJ , and use it to prove global convergence in certain classes of CRNs where 
only very mild assumptions are made on the kinetics. 

The key geometrical insight on which the results are built was provided in 
[B| and generalised in Q . Stated very briefly, these results show that sometimes 
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the existence of an integral of motion in a strongly order-preserving dynamical 
system allows the construction of a Liapunov function on each level set, which 
increases along nontrivial orbits. This in turn has strong implications for the 
convergence of orbits. While in general it may be an unusual conjunction of 
affairs to have a first integral, strong monotonicity, and moreover integral and 
order cone related in a particular way, it actually appears that this situation 
is not uncommon in CRNs. However identifying when this situation occurs 
is nontrivial, and explicit construction of families of CRNs satisfying all these 
conditions becomes important. 

The results at several points will be presented in considerably less generality 
than possible in order to simplify the presentation and highlight the key geo- 
metrical points, although where a theoretical result allows greater generality, 
this may be mentioned. 

2. A convergence result 

The result presented in this section, and applied in subsequent sections, is 
essentially derived from Theorem 2.4 in [6|], with slight modification and special- 
isation for our purposes. Note that Theorem 2.4 in Q stated that all orbits on 
a level set of the system in question converge to a unique equilibrium, provided 
the equilibrium exists, while what was actually proved was that all bounded or- 
bits converge to the unique equilibrium. As remarked in Q , it remains an open 
question whether the word "bounded" can be dropped from the statement of the 
main result as erroneously done in [6|. However for the purposes of this paper 
all orbits will be bounded and only convergence of bounded orbits is required. 
Note also that much of the difficulty in the proof of Theorem 2.4 in Q stemmed 
from the fact that the integrals of motion considered were in general nonlinear 
whereas here only the linear case is required. 

Standard notions from convex geometry (as in 0, 0, |T(| for example) will be 
assumed. Closed, convex and pointed cones in R™ define partial orders on R n . 
Following Q, cones in R n which are closed, convex, pointed and solid will be 
referred to as proper. Standard results in the theory of monotone dynamical 
systems [3, |4[ will also be assumed. 

Notation. The symbols <,>,<, ^, » will refer to the usual partial or- 
dering of vectors in R™ derived from the nonnegative orthant R> . So, given 
a,b e R n , a b means b — a e RJ , a < b means b — a e RJ \{0} and a « b 
means b — a e int(R> ). When the ordering is defined by some other closed, 
convex and pointed cone K, the alternative symbols <,>,<,>,-<,> will be 
used. So, for example, a -€ b means b — a e int K, and so forth. 

Let Y, K be proper cones in R™ with K => Y. Define K* to be the dual cone 
to K, i.e., K* = {y e R™ | (y, k) ^ for all k e K). Consider a system 

x = F(x) (1) 

on Y, and assume that ([1} defines a local semiflow <f> on Y. 
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2.1. Three conditions 

We define three conditions on (|T|) and the associated semiflow which will be 
referred to as Conditions 1, 2 and 3: 

1. tf> is monotone with respect to K, i.e., given x, y e Y with x < y and any 
t > such that <t>t{%) and <j>t(y) are defined, (j)t(x) < <t>t(y)- Moreover, <j) 
is strongly monotone in int Y in the following sense: given x < y with 
at least one of x or y in int Y, then 4>t(x) < <fit{y) for all t > such that 
0t(a;) and </>t(y) are defined. 

2. The system has an increasing linear first integral, namely there exists 
a linear function H:R n — > R such that (i) Vi? e int (!£"*), and (ii) for 
each y e Y, (VH,F(y)) = 0. Assume (w.l.o.g.) that H(0) = 0. 

3. 4> has no limit points on dY\{0}. 

2.2. Remarks on the conditions 

Condition 1. Verifying monotonicity/ strong monotonicity of 4> is generally 
the most nontrivial step in applications. The natural approach is via examina- 
tion of the Jacobian matrix of the system as will be spelled out in Section 12.41 
below. 

Condition 2. For each /ieK define the level set Sh = {y e Y \ H(y) = h). 
The condition (VH,F(y)} = tells us that H is constant along trajectories, 
i.e., the sets Sh are invariant. Choosing any set of basis vectors for (Vff ) x and 
writing these as the columns of a matrix T we can write F(y) = Tv(y) for some 
function v. Y — » R ,l_1 . The implication of Vif e mt(K*) is as follows: given 
x, y e Y with x > y, by convexity of Y , the line segment between x and y lies 
in Y, and by integrating VH along this line segment we get H{x) > H(y). So 
H is increasing with respect to the partial order defined by K. This implies in 
turn that the level sets of H are unordered, and since H(0) = 0, H{y) > for 
all y e Y\{0} (so So = {0}). Since H is linear and Y is a cone, clearly the range 
of H is [0, go), that is, H : Y — > [0, go) is surjective. 

Condition 3. Condition 3 is stronger than standard "persistence" assump- 
tions which require only that limit sets of interior points in Y may not intersect 
the boundary. The particular class of examples discussed here can be shown 
to fulfil Condition 3; however making a weaker persistence assumption leads to 
similar, slightly rephrased, convergence results to the ones presented here. The 
choice to make a strong persistence assumption is to reduce technical digressions 
and allow the main arguments to be more transparent. 

2.3. The key theoretical result 

The basic theoretical result from which applications to CRNs will follow is: 

Theorem 1. Consider System fl]) with Conditions 1, 2 and 3 above. For each 
h e [0, oo), Sh contains a unique equilibrium e(h), and e(h) attracts all of Sh- 

Proof. Let E = {y e Y F(y) = 0} be the equilibrium set of the system. The 
main landmarks in the proof are presented below, but the reader may refer to 
to fill in the details. 
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1. Any two equilibria ei, ei must satisfy e\ >- ei or e\ -< e%. It follows from 
Conditions 2 and 3 that E n <5Y = {0}. Consequently given two distinct 
equilibria ei,e2 at least one of these must lie in intF. As a result, if 
e~i,ei satisfy ei < e2, then by strong monotonicity in intF, ei -< ei. Now 
suppose there are equilibria e±, e~i which are unordered by <. By repeating 
the geometrical arguments as in the proof of Lemma 5.11 of [g we obtain 
the existence of an equilibrium e such that e < e~\ (say) but not e -€ ei, 
contradicting the fact that e < e\ implies e -< e±. 

2. Each Sh is a compact, convex set: by linearity of H, Sh is the intersection 
of a hyperplane with Y and so convex, while boundedness follows from 
the condition that Vif e int(iT*) (see Lemma 3.9 in Q). Thus each level 
set contains an equilibrium. Uniqueness of the equilibrium follows from 
the fact that Sh is unordered, whereas from 1, two equilibria e\ and ei 
must satisfy e\ > ei or e\ -< e%. 

3. The map e : [0, oo) — > E which associates to each h the unique equilibrium 
e(h) on Sh is an order-preserving homeomorphism. That e is an order- 
preserving bijection with continuous inverse (namely H\ E ) is immediate 
from parts 1 and 2, and linearity of H. That e is continuous follows 
essentially from the fact that E is closed (see Lemma 5.12 in Q). 

4. There is a continuous scalar function L : Y — > [0, oo), increasing along non- 
trivial orbits in int Y, and which attains a global maximum on Sh at e(h). 
The proofs of these facts follow the related results in (Lemmas 5.13 
to 5.18). The value of L(y) is the value of H at the unique point Q{y) 
where y — dK intersects E. Lemmas 5.13 to 5.17 in Q carry over to the 
present case immediately. We claim that L is increasing along orbits of 
those y $ E for which L(y) > 0. Indeed, consider any such y. By the def- 
inition of Q(y), Q(y) < y. Since L(y) = H(Q(y)) > 0, we have Q(y) ^ 0, 
so, by Condition 3, Q(y) e intF. Strong monotonicity in intY implies 
that for any t > 0, <f>t(Q(y)) = Q(y) < <M?/)- Consider any e e E. If 
e < Q(y) then e < <f> t (y), i.e., e # Q(Mv))- So Q(Mv)) > Q(v), and 
thus L(<p t (y)) > L(y), which proves the claim. It remains to notice that 
L(y) > for any y e intl", since otherwise y e which contradicts 
int y c int K. We remark that L is well defined at points on dY, but may 
not be (strictly) increasing if L(y) = (this is possible if y e dY n dK) as 
the line segment [0, y] lies in where we do not necessarily have strong 
monotonicity. 

5. If h = it is trivial that e(h) attracts all of S(h) as Sh = e(h) = {0}. 
So assume h =(= 0. Since 5^ is compact, each x e Sh has a nonempty w- 
limit set, ui(x). The Liapunov function in intK, strictly increasing along 
nontrivial trajectories and achieving a maximum at e(h), guarantees that 
lj(x) does not intersect int Y\{e(h)}. So lj(x) £ BY u {e(/i)}. But by 
assumption cj(a;) does not intersect <9F\{0} and ^ 5/i, so a; (a;) = {e(h)}. 

□ 

Remark. For h =(= the equilibrium e(/i) defined in the above proof clearly 
lies in inty and so in the relative interior of Sh- 
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2.4- Jacobian conditions for monotonicity and strong monotonicity 

In general, the greatest difficulty in applying Theorem [T] is verifying mono- 
tonicity of a local semiflow in Y and strong monotonicity in int Y in the sense 
described in Condition 1. For this reason it is useful to have conditions which 
can more easily be checked for a system. 

Terminology. Let J be a square matrix. Given a closed, convex and pointed 
cone K E R™, we will say that J is JC-positive if y e K implies Jy e K. A K- 
positive matrix J is X-irreducible if the only faces F of K for which J{F) E F 
are {0} and K . J is ^-quasipositive if given any y e K, there exists aei 
such that Jy + ay e K. Suppose that there exists a e R such that J + al is 
-FT-positive and ^-irreducible, then we say that J is strictly i-T-quasipositive. 

The following is a slight adaptation of results in Q for our purposes. 

Lemma 2. Consider a proper cone K E R™ ; some open set U E R™ , and 
a C 1 vector field f:U—> R™ with Jacobian matrix Df. Let X c U be some 
convex domain, positively invariant under the local flow 4>u defined by f , and 
let (f> be the induced local semiflow on X. Assume that Df is K -quasipositive 
in X. Consider some Xq,Xi e X with xq < X\. Then 4>t(xo) < 4>t(xi) for each 
t > such that 4> t (xo),(j)t(xi) exist. If there exists yo on the line segment joining 
xo to X\ such that Df(yo) is strictly K -quasipositive, then <j>t{xo) -< 4>t{x\) for 
each t > such that 4>t{xo),4>t(xi) exist. 

Proof. By basic results on monotone dynamical systems (see Section 3.1 in 
[4]), if-quasipositivity implies monotonicity: given Xq,X\ e X with xq < X\, 
4>t(xo) < ^(^l) f° r each t > such that (f>t(xo),4>t(xi) exist. Note that if 
Df(yo) is strictly -FT-quasipositive, then for each x e dK \{0}, there exists v 6 K* 
such that v(x) = and v(Df(y )x) > (see Proposition 3.10 of \4f). The final 
statement now follows from Lemma 3.7 in □ 

Remark. Since we are discussing the Jacobian matrix it reduces technical- 
ities to assume that the vector field is defined on some open set containing the 
positively invariant set X in which we are interested. 



3. Some background material 

From now on we fix the state space of interest as R™ > the nonnegative 
orthant in R m and focus on systems with a linear first integral. Let T be 
an to x n matrix, and v : U — > 1" be a C 1 function where U is some open 
neighbourhood of R™ - We consider the system 

x = Tv(x) (2) 

defining a local flow </> on U. We assume that R™ is (positively) invariant under 
<j> and so get a local semiflow tj> on R Jo- 
Dynamical systems arising from chemical reaction networks typically take 
the form In this context, the matrix V is termed the stoichiometric 
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matrix of the system, and the intersection of any coset of Im T with E™ is called 
a stoichiometry class. Clearly stoichiometry classes are positively invariant 
under <j>. The function v describes the rates of reaction. A variety of common 
forms for v exist in the experimental and modelling literature, all of which satisfy 
the weak assumptions we present later fSection l5.ip . 

3.1. Qualitative classes and sets of matrices 

A real matrix M determines the following three sets of matrices: 

1. The qualitative class Q(M) [ljj of all matrices with the same sign pattern 
as M. Explicitly, N e Q(M) if N has the same dimensions as M and 
satisfies sign(Ay) = sign(My), i.e., My > => Ay > 0, My < => 
Nij < and My = Q =^ Ay =0. 

2. The set Q Q (M) = cl(Q(M)), the closure of Q(M). Q (M) consists of 
all matrices A such that My > =*• Ay 3j 0, My < => Ay s$ and 
My > % = 0. 

3. The set Qi(M) consisting of all matrices A such that My > => Ay ^ 0, 
M„- • : A^y • 0. 

The following fact will often be applied without comment: if c e Qo(b) and 
b e Qi(a), then c e Qi(a). The next lemma can easily be checked directly. 

Lemma 3. Let a, b e R™, and Zei D be an n x n diagonal matrix with positive 
diagonal entries. Suppose Da,Db are (—1,0, l)-vectors. Then a + b e Qi(a) n 
Ci(6). □ 

5.1!. Associating digraphs with CRNs 

A DSR graph ("Hi is a signed, labelled, bipartite, multidigraph, examination 
of which can lead to a number of conclusions about the behaviour of a system 
of differential equations. In this paper we will be interested only in whether a 
given digraph is strongly connected, namely whether there is a (directed) path 
from each vertex in the digraph to each other vertex. So, following [13|, we 
present a closely related but simplified construction. Given an m x n matrix 
A and an n x m matrix B, we construct a bipartite digraph Ga,b on m + n 
vertices as follows: associate a set of m vertices Ui, . . . , u rn with the rows of A; 
associate another n vertices Vi, . . . ,v n with the columns of A; add the arc UiVj 
iff Aij =|= 0; add the arc VjUi iff Bji =(= 0. 

Such a digraph G(x) can be associated with System ^ at each point in 
x e M™ as follows: define V(x) = [<5«j(a:)/da;j], and let G(x) = C?r,v(a:)- 

Remark. It is in general possible to factorise a given function in different 
ways, say F(x) = Tv(x) = Tv(x). These different factorisations result in dif- 
ferent digraphs. Thus the digraph associated with System ([2]) is not strictly 
associated with the vector field itself, but with a particular factorisation see 



12j for more discussion of this. In the case of CRNs, there is a natural factori- 



sation, where T is the stoichiometric matrix. 
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3.3. The n-cube 

Given a matrix M, Mk will refer to the fcth column of M, and M k to the 
fcth row of M. 



Fix n > 2 and define the nx2™ matrix B as the (0, l)-matrix whose jth 
column read vertically upwards as a binary string gives the number j ' — 1 (so 
2r=i 2 .By — 3 ~ !)• For example, if n = 3: 

10 10 10 1 
Zi =100110011 
1 1 1 1 

Note that the columns of B are just the vertices of the n-dimensional hypercube, 
and two vertices are adjacent if they differ in exactly one coordinate. The choice 
of ordering is important: given i < j, Bi and Bj are adjacent iff (i) there exists 
k > such that j = i + 2 k and (ii) i - 1 < 2 k mod 2 k+1 . The following is 
immediate: 



, ifi-1^ 2 k - 1 mod 2 k , 
= e k . □ 

The geometrical meaning of Lemma 2] is illustrated for n = 3 in Figure [1] 



Lemma 4. Given i e {1, . . . , 2™} and k e {1, 

i/ien Z?i — B i _ 2 k-i = efc cmiZ otherwise B i+2 k- 



Bi 



X3 



/ 



X-2 



/ 



/ 



/ 



- 4 
/ 

2 >xi 



Figure 1: The 3-cube in R 3 with vertex n corresponding to column not B. Bkn is the (4 — fc)th 
entry in the binary representation of n — 1. Observe that if i — 1 < 2 fe_1 mod 2 , then 
B j+2 k-i — Bi = corresponding to the geometrical fact that each basis vector is associated 
with 4 parallel edges of the cube defined by B. For example B 3+2 2 ~ E>3 ~ ^7 ~ ^3 — 
(0, 1, 1) T - (0, 1, 0) T = (0, 0, 1) T = e 3 . 



3.4- Cubic cones 

We now construct a class of cones which we refer to as "cubic cones" . First 
we show that no binary vector is in the convex hull of other, different, binary 
vectors. Geometrically this is obvious: these binary vectors are precisely the 
vertices of the n-cube. 

Lemma 5. Let B be any (0, 1) -matrix without repeated columns and let p be a 
nonnegative vector with Y^iPi = 1- If Bp = q and q is a (0, V)-vector thenp = ej 
for some j. 
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Proof. We look for solutions to q = Bp with q a binary vector and XiiPi = !• 
For each i, let Sf be the set of indices such that k e Sf Bi k = 0. Similarly 
define S} via k e S} o B lk = 1. Then: 

1. If q { = Xifc -BifcPfc = 1, then J] kES i Pk = 1, so 2fc 6 g9 Pfc = 0, so p fc = for 
each fc e 5°. 

2. If qi = 2 fc SifePfc = 0, then 2 fceS i Pfc = 0, so p k = for each ke S}. 

Since p is not the zero vector, there exists j such that pj =|= 0. Then for each i 
such that qi = 1, j $ Sf, i.e., j e , i.e., -By = 1; similarly for each i such that 
<7i = 0, j $ Sf, i.e., j e S?, i.e., B,j = 0. We have g< = By, i.e., g is just the 
jth column of B. If there exists I =(= J such that pi =(= 0, then we get <7i = _Bi;. 
Since B has no repeated columns this is impossible, so pi = for i =(= j- Since 
SiPi = 1) P = e i- □ 

The construction in the next lemma is crucial to what is to follow: 

Lemma 6. Let m > n and consider an m x n matrix T of rank n. Let c e K m 
be some vector not in ImT. Define 1 to be the vector in M 2 each of whose 
entries is 1 and define A = cl T + TB. With A defined in this way, 

(i) K = {Az | z > 0} is a closed, convex and pointed cone in R m ; 

(ii) The set of those r e kerT T for which r T c > is nonempty and is contained 
in int(isT*); 

(hi) The nonnegative multiples of the 2™ columns of A are precisely the closed 
one- dimensional faces of K . 

Proof. That K is a closed, convex cone is immediate from the definition. Since 
c ^ Iml, there exists r e kerT T such that r T c > 0. For any p > and any such 
r, r T Ap = r T cl T p = r T c^pi > 0. Sore int(K*). Since int(i4T*) is nonempty, 
K is pointed. 

The third statement fails to be true only if some column of A can be written 
as a nonnegative combination of other columns of A. Suppose Ae& = Ap where 
0<p. Then: 

(cl T + TB)e k = (cl T + TB)p 

i.e., 

c(l-Y J P l )=TB(p-e k ) 

Left-multiplying both sides by r T gives r T c(l — = 0> so YiPi = 1 an d 

TB(p — e k ) = 0. Since T has rank n, kerT = {0}, and so B(p — e k ) = 0. By 
Lemma O p = e k . So each column of A is an extremal vector of K. □ 

Geometrical interpretation. Since T defines a linear map from R™ to 
R m with trivial kernel and B has rank n, the columns of TB are the extremal 
points on an n-dimensional parallelotope in R m . Adding on the constant vector c 
simply translates these points, so the columns of A also define an n-parallelotope 
in K m . Ensuring that the constant vector c does not lie in ImT means that the 
cone K{A) = {Az \ z > 0} is pointed, and the convex hull of columns of A, 
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conv(A), is a cross-section of this cone. As conv(A) is simply the n-cube after 
an affine mapping we call these cones cubic cones and the associated partial 
orders cubic orders. 

Obviously, the incidence structure of conv(A) is precisely that of the cube, 
conv(B), namely the convex hull of any 2 k vectors A ri , . . . , A r ^ k is a fc-dimensional 
face of conv(A) if and only if the convex hull of the corresponding vectors 
B ri , . . . , B r2k is a fc-dimensional face of the n-cube. 

Each vector is associated with 2 n_1 edges of conv(A) corresponding to 
2™~ 1 parallel edges of conv(-B). More precisely, by Lemma [H if i — 1 ^ 2 k ~ 1 
mod 2 k then 

A, - A t _ 2 k-i = c + TBi - (c + rs 2 _ 2fc -i) = r(Bi - Bi_ 2k -i) = T k 
while if i — 1 < 2 k ~ 1 mod 2 k then 

A t+2k -i -Ai = c + TBi - (c + TBi_ 2h -i) = T(B l+2 ,-i - B t ) = T k 

So for any k e {1, . . . , n] and any i e {1, . . . , 2"} we define 

... . = J i - 2 & - 1 (if i — 1 5s 2 fc - 1 mod 2 k ) 
J(h > ~ { i + 2 k ~ 1 (otherwise) 

and can always write for any i, k 

T k = sign(j(i, fc) - i)(A j( ^ k) - A t ). 

4. Systems preserving cubic orders 

We proceed to construct a class of systems to which Theorem [1] can be 
applied. We do this by prescribing a priori the geometric nature of the order 
to be preserved and then constructing systems which preserve such orders. Via 
the next two theorems we provide a set of easily checkable conditions which 
ensure that System ^ has the behaviour described in Theorem [T] Later we 
will provide examples of CRNs which fulfil these conditions. 

Theorem 7. Consider System (0) with T now an (n + 1) x n matrix of rank n. 
Assume: 

1. For some c =(= 0, A = cl T + TB satisfies: 

(a) There exists a nonnegative right-inverse to A, i.e., a nonnegative 
matrix P such that AP = I . 

(b) There exists a diagonal matrix D with positive diagonal entries such 
that DA is a (—1, 0, T)-matrix 

2. For each x e V(x) = [dv l {x)/dx J ] e Q (-T T ) 

3. For each x e im^R^Q 1 ), the digraph Gry( x ), is strongly connected. 

Define 

K = K(A) = {Az\z^ {)}. 

Then: 
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1. MJo 1 S K. 

2. c ^ ImT. in^Jf*) is nonempty and K is pointed. 

3. Each column of A is an extremal vector of K. 

4. For each x e R™^ 1 , TV(x) is K-quasipositive, while for eachx e int(R>g ), 
TV{x) is strictly K-quasipositive. 

Before proving Theorem [7] we show how when it is combined with Theorem[T] 
we get a global convergence result. 

Theorem 8. Consider System (0) with the assumptions in Theorem^ Assume 
further that 5RJq 1 \{0} contains no limit points of the local semiflow <p. Then 
there exists a linear function H : KJq 1 — » R whose level sets 

S h = {ye RJ+ 1 | H(y) = h). 

are invariant under <f>. Moreover if for some h, Sh is nonempty, it contains a 
unique equilibrium e{h) which attracts all of Sh- 

Proof. The cone K in Theorem [Jj is clearly closed and convex. That it is 
pointed and solid follow from results 1 and 2 of the theorem. So if is a proper 
cone. Result 4 of Theorem [JJ gives that the Jacobian matrix of the system is 
if-quasipositive everywhere in RJq 1 and strictly if-quasipositive in int(R>Q ). 
By Lemma[21 if-quasipositivity implies that given xq <x%, (f>t(xo) < 4>t{x{) (for 
each t > such that (/>t(xo),<frt(zi) exist); further, strict i\-quasipositivity in 
im^R^Q 1 ) implies that if either of xo or x\ lies in miiWt^ 1 ) then <j>t(xo) -€ 4>t{xi). 
Thus Condition 1 of Theorem Q] is satisfied. 

Since by Theorem[7j2), c ^ ImT, LemmaEJii) allows us to choose r e kerr T 
such that r T c > 0. Define H(x) = r T x. H is by definition linear, with H (0) = 0, 
and clearly 

(VH,Tv(y)) = r T Tv(y) =0 

so that level sets of H are invariant. Lemma [6jii) provides further that VH = 
r e int(K*), so the integral is increasing. Thus Condition 2 is satisfied. 

Finally, Condition 3 is assumed. Thus the assumptions of Theorem [T] are 
fulfilled and so Theorem [1] can be applied. □ 

A somewhat more general statement than the following is proved as Theo- 



rem 1 in 13 1 



Lemma 9. Let K E R" be a closed, convex and pointed cone, A an n x m 
matrix, and B' anmxn matrix. Suppose that Im A $ span F for any nontrivial 
face F of K, and that AB is K-quasipositive for each B £ Qq(B / ). Then 
whenever G a,b is strongly connected, AB is strictly K-quasipositive. 
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4 ■ 1 ■ Proof of theorem [7] 

The lemmas presented so far make the proof now straightforward. 

Proof of Theorem El 1. M^ 1 c A. Since APei = e, and Pei is nonnega- 
tive, ei e K for each i and so by convexity MJq 1 £ K. 

2. c ^ ImF. There exists r e int(A*) and A is pointed. From the 
previous part, R n+1 = K — K. One consequence is that, regarded as a map 
from R 2 to R n+1 , A is surjective. It follows that c ^ ImT: otherwise, since all 
vectors in Im A are of the form Ac + Tb (for some scalar A and some vector 6), 
each vector in Im A could be written Fd for some vector d; however T only has 
rank n contradicting the fact that A is surjective. 

Since c ^ ImF, by Lemma [BJii) there exists r e kerT T for which r T c > 0, 
and r e int(AT*). Thus if is pointed. Note that since MJq 1 c A", it follows that 

int(A*) c int(M^o X )i so r » °- 

3. Each column of A is an extremal vector of A. This follows from 
Lemma IHIiii) • 

4a. For each y e K™^ 1 , TV{y) is A-quasipositive. Fix y and let V = 

V(y). Recall that for i e {1, . . . , 2™} and k e {1, . . . , n} there exists k) e 
{1, ... ,2"} such that 

T k = sign(j(«, k) - i)(Aj(»,fc) - AO- 
Define = T^ fc Ai|, cti = 2 fc ait and a = max^ ctj. We have 

k 

If i > j(i,k), then = Aj — A^m, and by Lemma H we have e Qx(Aj). 
Since -F fe £ Qo(r^), -V fc e Qi(Af), so -V k A t > 0. Similarly, i < j(i,k) 
implies V k Ai ^ 0. So, in either case 

T k V k A l + a lk A, = sign(j(i, k) - i)(A i(a) - Ai)V k Ai + a ik Ai 
= (Aj(i,fe) - Aj)a ife + a 4 fcAi = Oik^-j(i,k)- 

So 

+ a 4 A, = '£(r k V k A i + a lk A t ) = £ a ife A i(i>fe) £ AT (3) 
fe fc 

since ^ for each i. Given any x = ijAj e A (ti > for each i): 
TVx + ax = ( aA ' ; + rvA >) 

i 

= J^Uda-a^Ai+TVAi + aiAi) (4) 

i 

Since a — on > for each i and YV Ai + c^A,; E A, TVx + ax e A. 

4b. For each y e int(RJ ), ry(y) is strictly A-quasipositive. By 

4(a) we already know that TV is A-quasipositive for all V in Qo(~ r T ). By 



11 



assumption, for each y e int(M^Q 1 ), the digraph G r y( y ), is strongly connected. 
In order to apply Lemma [9] we need to confirm that Im T span F for any 
nontrivial face F of K. 

Consider a nontrivial face F of K. We put the columns Aj e F into a matrix 
Ajr. Suppose ImT £ ImAi?, i.e., there exists a matrix T such that T = ApT. 
Since A^ simply consists of a subset of the columns of A, for any i there exists 
such that Apei = c + TBjrq, and as T = ApT, c = Ap(e-i — TBju\). Thus 
A = A f [(e, - TBj^)l T + TB]. This shows that ImA £ ImA F , contradicting 
the fact that F is a nontrivial face of K. Thus ImL $ spanF for any nontrivial 
face F of if. 

It now follows from Lemma[9]that TV(y) is strictly if-quasipositive for each 
yeint(R5 ). ' □ 



5. Applications to chemical reaction networks 

5.1. Assumptions on reaction kinetics: Condition K 

We consider chemical reaction networks of the form defined on M™ . At 
various points above we have mentioned weak assumptions on the kinetics of 
reactions, and here we spell these out. Let xi refer to the concentration of the 
ith chemical. Each reaction has a left-hand side and a right-hand side: let Xjj 
be the indices of chemicals occurring on the left of reaction j and Z,. r be the 
indices of the chemicals occurring on the right of reaction j. Define L^ to be 
the difference between the number of molecules of the zth chemical occurring on 
the right of reaction j and on the left of reaction j. Let Vj represent the rate at 
which reaction j proceeds to the right. With these definitions we indeed get the 
dynamical system ([2|) for the evolution of the concentrations. We assume that 
reactions are either reversible or irreversible, and for an irreversible reaction we 
make the arbitrary choice that the substrates occur on the left and products on 
the right so the reaction always proceeds with nonnegative velocity to the right. 

We assume: 

(Kl) Tij(dvj/dxi) sS 0, and if L^- = then cvj/dxi = 0. In other words 

V(x) = dV g^ e Qo{— r T ). This condition has been discussed before in 

[lH HH for example, and is satisfied by all reasonable kinetics provided no 
chemical occurs on both sides of a reaction. 
(K2) If reaction j is irreversible then: 

(i) Vj ^ with Vj = if and only if xi = for some i e 

(ii) If Xi > for all i e Xjj, then dvj/dxi > for each i e Ijj. 
(K3) If reaction j is reversible then: 

(i) If Xi = for some i e Xj i (resp. for some i e 2}, r ) then Vj < (resp. 
Vj >0). 

(ii) If Xi = for some i e Xj t i (resp. for some i e Ij, r ), then Vj < (resp. 
Vj > 0) if and only if x# > for each i' e Xj r (resp. for each i' e Xj,i). 

(hi) If Xi > for all i e Xjj (resp. for all i e Tj r ) then, for each i e Ijj 
(resp. i elj, r ), dvj(x)/dxi > (resp. dvj{x)/dxi < 0). 
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We will refer to these very mild and natural assumptions collectively as 
Condition K. The following simple result will be useful in the sequel: 

Lemma 10. Let a system of the form 0] satisfy Condition K. Then for any x e 
R™ > an y 3 an d an y i such that Xi = there holds TijVj(x) > 0. Consequently, 
for such a system R™ is positively invariant. 

Proof. Suppose the contrary: then there exists a reaction j, some x e R™ 
and a chemical i such that Xi = and TijVjix) < 0. First of all, if chemical i 
does not participate in reaction j then = 0, so this cannot be the case. If 
i occurs on the left side of j and j is irreversible then, by (K2)(i), Vj(x) = 0, 
consequently TijVj(x) =0. If z occurs only on the right of irreversible j then 
Tij > 0, and, again by (K2)(i), TijVj(x) ^ 0. Assume now that reaction j is 
reversible. Since, by application of (K3)(i), if chemical i occurs on both sides of 
reaction j then Xi = implies Vj(x) = 0, this is excluded. Consequently, either 
(i) chemical i occurs only on the left of reaction j and Vj (x) > 0, or (ii) chemical 
i occurs only on the right of reaction j and Vj(x) < 0. These possibilities are 
ruled out by assumption (K3)(i). □ 

Given a reversible system of reactions satisfying Condition K, knowledge of 
the sign pattern of T allows us to determine the sign pattern of V{x) at each 
x e R™ ; this in turn allows construction of the digraph at each x. (When some 
reactions are irreversible, if we have knowledge of how products affect reaction 
rates, then again we can determine the zero pattern of V(x) from T alone, and 
hence construct the digraph.) In particular, if all reactions are reversible then 
in fact for each x e int(R™ ), V(x) e Q(— r T ), and G TtV i x \ = G r>T T. 

5.2. Repelling faces of 'M> 

Let Sc {1, . . . , to} be a subset of {1, ... , to}, and S c = {1, . . . , m}\S. Define 

F s = {x e R m | x t > 0, i e S and x t = 0, i e S°}. 



Following [16|, Fs will be referred to as an elementary face of R™ . (Elemen- 
tary faces are the relative interiors of the closed faces of R™o-) If 5 is a proper, 
nonempty subset of {1, ... , m}, then Fs will be termed a nontrivial elementary 
face of R™ - Given a vector field on R™ defining a local semiflow on R™ i Fs 1S 
repelling if, at each x e Fs there exists i e S° such that Xi > 0. The following 
is almost obvious: 

Lemma 11. Consider a vector field on R™ which defines a local semiflow 
o?i R> , and assume that some elementary face F £ 5R™ is repelling. Then F 
contains no uj -limit points of (j). 

Proof. By the definition of a repelling face, each point on F is a start point 
of cf>: given x e F since there exists i such that Xi = and Xi > 0, there can 
clearly be no negative continuation of the trajectory through x. But since R^?o 
is locally compact , no w-limit set of <f> can contain a start point of <j> (Theorem 



5.6 in [17[) and so F cannot intersect any w-limit set of 4>. □ 
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Remark. The following implication of Condition K will be needed below. 
Given an elementary face Fs, some y, y' e Fs, some i e S c , and any j such that 
Tij =|= 0, then Vj(y) > (resp. Vj(y) < 0, resp. Vj(y) = 0) implies Vj(y') > 
(resp. Vj(y') < 0, resp. Vj(y') = 0). 

Lemma 12. Consider a system of the form (0j defined on ]R™ and satisfying 
Condition K. If for some elementary face Fs £ R™ there exists x e Fs where 
f(x) is not tangent to Fs (i.e., f(x) ^ sp&nFs) then Fs is repelling. 

Proof. If f(x) is not tangent to Fs, then there must exist i e S c with £4 ^ 0. 
As 

x i — 2} ^ij v j( x ) ano - eacn °f the summands is, by LemmaJTUl nonnegative, 
it follows that Xi > 0. As a consequence of the remark above, for any y e F s , 
Vi = Sj r ij v j(v) > °- Thus F s is repelling. □ 

5.3. A family of globally convergent CRNs 

We consider the following class of CRNs with n reactions onn+1 chemicals 
(n 5= 2). We refer to the nth member of this class as IZ^ and its stoichiometric 
matrix as T^h The CRNs are 



Fixing n, and choosing a convenient ordering for the reactions the stoichio- 
metric matrix V = T^ n ' is defined by the rules: 

1. Column 1 



ft< 3 > : 



A^B + D 
B^C^D 



A^B + E 
B ^ D ^ E 



(5) 



j = 1 




i = n + 1 
otherwise 



2. Columns 2, . . . , n — 1 



J = 2, • . • , 



n-1 




1 



i = n — j + 1 
i = n — j + 2 
otherwise 



3. Column n 




i = 1 
i = 2 
j = n + 1 
otherwise 



j = n: < 



1 i 
1 i 



r(«) for n = 2, 3, 4 are 




/ 1 \ 
1 -1 
-1 -1 
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-(4) 



V 





1 

1 -1 

-1 





1 \ 

-1 




1 ) 



Theorem 13. Assuming that all reactions are reversible, and the kinetics ful- 
fils Condition K, the CRN TZ^ (k > 2) has the following behaviour: each 
stoichiometry class contains a unique equilibrium to which all trajectories are 
attracted. 



Proof. The result follows from Theorem [8] if all the conditions in Theorem [7] 
can be shown to be met, and additionally, 5M>g \{0} contains no limit points 
of the local semiflow. Note first that for each n, T^ n > is an (n + 1) x n matrix 
and it is straightforward that kerl^™) is trivial. We fix n ^ 2, and for notational 
simplicity omit the superscript (n) from all objects. Defining B as before, let A 
be the (n + 1) x 2™ matrix: 



A 



l l T +TB. 



For example, for n = 3, 



A 



/ \ (1,1,1,1,1,1,1,1) 

1 

V o / 



I 



/ 



1 

V o i 



1 

1 -1 
-1 1 

1 -1 



+ 





1 

-1 



1 

-1 






1. Note that {Aj} are of the form 




ieS 



10 10 10 1 
1 1 1 1 
1 1 1 1 



where S can be any subset of {1, . 



We have: 



e n = Ai, 
e n+ i = e 

&n — k 

ei = e n + e 



Ti = A 2 



A 2 fc + i, k 



1. 



f=2 i , — + n, — i, . . . ,n • 

i + S?=2 r i = A 2 + A 2 „_ 1 . 
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Each of these facts can quickly be proved by induction. Consequently, 
there exists a nonnegative matrix P such that AP = I. For example, for 
n = 3, 



P = 



f° 





1 




1 








1 





1 












































1 











V« 











Thus Condition 1(a) of Theorem [7] is met. 

2. It is easy to confirm that the kth entry of any sum of columns of T has 
value —1, or 1 for k =j= n and value —2, —1 or for k = n. Consequently, 
A is a (—1,0, l)-matrix. Thus Condition 1(b) of Theorem[7]is met. 

3. An immediate consequence of Condition K is that for each x e M™^, 1 , 
V(x) = [dvi{x)/dxj] e Q (-T T ). Thus Condition 2 of Theorem [7] is met. 

4. Using (K3)(iii), the digraphs G-py{ x ) a t any point x e in^R™^ 1 ) can be 
inferred. These are illustrated for n = 2,3,4 in Figure [51 Clearly they 
are strongly connected, each consisting of a single double-cycle with one 
offshoot. Thus Condition 3 of Theorem [7] is met. 

5. Each nontrivial elementary face of KJq 1 * s repelling. If this is not so for 
some nontrivial elementary face Fs £ IR^q 1 , then by Lemma fT2l Tv(x) is 
tangent to Fg at each x e Fg. This in turn implies for each x e Fs and 
each i e S°, that i< = 0. Since, by Lemma [TOl i e S c implies TijVj(x) 5= 
for each j, Xi = now implies that TijVj (x) = for each j. Thus whenever 
r,-j =|= then Vj(x) = 0. Consider some j such that =(= 0. Note that no 
chemical occurs on more than one side of this reaction system and that 
all reactions are reversible. Suppose that chemical i appears only on the 
right of reaction j so that IYj > 0. Vj(x) = then implies, via (K3)(ii), 
that there exists i' e Ijj n S c , and so T^j < 0. Suppose, on the other hand 
that chemical i appears only on the left of reaction j so that Fy < 0. In 
this case, via (K3)(ii), there exists i' e Xj >r n S c and so T^j > 0. Thus, for 
a system of reversible reactions such as this, with no chemical appearing 
on more than one side of any reaction, the condition that Tv(x) is tangent 
to Fg at each x e Fs has a simple interpretation in terms of T: each 
column of the submatrix T s of T with only rows indexed by S c must 
contain either both a positive and a negative entry, or only zeros. That 
this does not occur in r*™) for any nonempty S c {1, . . . ,n + 1} can be 
checked. If S includes 1, it must include 2 or n + 1; if <S* includes k where 
2 ^ k < n, then it must include k + 1 and k — 1; if S includes n + 1, 
then it must include 1 : thus if 5* is nonempty then S = {l,...,n + l}. 
In the terminology of [18|, the systems have no "siphons". In physical 
terms, setting any proper subset of reactant concentrations to zero, there 
always remains some active reaction producing one of the absent reactants 
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with positive rate. Consequently, by Lemma [TT1 no nontrivial elementary 
face of MJq 1 contains any w-limit points of the associated semifiow <f>. As 
3R>g \{0} is precisely the union of these nontrivial elementary faces, it 
contains no w-limit points of the semifiow. 



□ 




E n = 4 C 



A . ^N, 



I " = 3 1 

D C \ X 

Figure 2: The digraphs associated with reaction networks in l(5} for n = 2,3,4 under the 
assumption that ^ ^^f~ j e S(— T T )- Vertices associated with rows of T (i.e., with particular 
chemicals) are labelled, while those associated with columns (i.e., with reactions) are filled 
circles. Antiparallel edge-pairs are shown as bidirected edges. As the only matter of interest 
is connectivity of the graphs, edge signs and labels are omitted. 



5.4- Another example 

The family of CRNs presented in the previous section are not the only CRNs 
preserving cubic orders. As another example of a CRN which preserves a cubic 
order, consider the following: 



2A ^ B ^ C + D 
C^D 



(6) 



giving 

f 2 \ 

1 -1 

-1 -1 

V i -i 0/ 

kerT = {0} and the positive vector (1, 2, 1, 1) T lies in kerT T . Defining B in the 
usual way, 



A^esl 1 +TB 



( 



1 

Vo i 





1 



-1 



2 2 

-1 -1 

1 

1 



2 \ 


-1 
o / 
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Left-multiplying A by the diagonal matrix diag{l/2, 1, 1, 1} gives a (—1,0,1) 
matrix. The matrix 



/ 



1/2 







1/2 



\ 



can be confirmed to be a nonnegative right-inverse for A. Consequently, assump- 
tion 1 in Theorem [7] is satisfied. Assuming that the kinetics fulfils (K3)(iii), the 
digraph in int(R|. ) is shown in Figure [3] and is clearly strongly connected. 



• C 

I 3 



Figure 3: The digraph Gr.v for the reaction system in Conventions are as in Figure [2] 

Thus all the conditions of Theorem [7] are fulfilled. Further, it can easily be 
checked, using the approach in part 5 of the proof of Theorem Q21 that each 
nontrivial elementary face of MJg 1 is repelling and so by Lemma HT1 5RJg 1 \{0} 
contains no w-limit points of the associated semiflow <f>. Thus by Theorem [8] 
each stoichiometry class contains a unique equilibrium to which all trajectories 
are attracted. Clearly with the exception of the equilibrium at 0, these globally 
attracting equilibria lie in vatiWi^ 1 ). 

We remark that in the above examples, the assumption that all reactions 
are reversible can be relaxed provided that this is done in such a way that 
the digraph in int(MJ+ 1 ) remains strongly connected, and such that nontrivial 
elementary faces of KJg 1 remain repelling. 

6. Concluding remarks 

Systems of differential equations arising in chemistry and displaying strong 
convergence properties have been presented. The examples are only intended 
to illustrate the applicability of the main theorems (Theorems [7J and |5]) : char- 
acterising fully the set of CRNs for which global convergence can be proved via 
these theorems remains a task for the future. 

In addition, a number of variants on the same theme remain to be explored. 
While this work has focussed on systems preserving so-called cubic orders, these 
are by no means the only partial orders which may be preserved by a system 
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of chemical reactions as demonstrated in [19]. Closely related results can be 
derived where the geometrical construction in this paper is restricted to linear 
invariant subspaces: this approach throws up further examples of chemical re- 
action networks with strong global convergence properties. These themes will 
be explored in forthcoming papers. 
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